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Abstract. We review recent results on properties of the meson gas relevant for Heavy Ion Collision 
and Nuclear Matter experiments, within the framework of chiral lagrangians. In particular, we 
describe the temperature and density evolution of the a and p poles and its connection with chiral 
symmetry restoration, as well as the chemical nonequilibrated phase and transport coefficients. 
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HOT AND DENSE RESONANCES IN UNITARIZED CHIRAL 
PERTURBATION THEORY 

It is well established nowadays that QCD undergoes a chiral restoration transition, as 
confirmed by different lattice simulations [1,2]. The transition is presumably a crossover 
for 2+1 massive light flavors, which means that the critical point could be different for 
different observables. In the lattice, the main order parameters are the quark condensate 
and the light scalar susceptibility. Nevertheless, there have been several efforts in the 
literature to try to find signals of chiral symmetry restorations in physical observables 
which could be eventually measured in Heavy Ion or Nuclear Matter experiments. 
The most promising signals are related to the behaviour of resonances in hot and 
dense matter, such as the p — a\ degeneracy [3, 4] or the role of the p resonance in 
the dilepton spectrum [3, 5]. In the latter case, a significant dropping of the p mass 
would be interpreted as a signal of chiral restoration according to the Brown-Rho 
scaling [6] or QCD sum rules [7] scenarios. However, this does not seem to be the case 
experimentally in Heavy Ion Collisions, where the alternative explanation of broadening 
dominance [3] is equally valid to explain the data and it is actually favored in recent 
experiments, such as NA60 with dimuons [8]. At RHIC energies (PHENIX), even the 
broadening picture is insufficient to explain the observed excess in the region of low 
invariant mass [9, 10]. The situation is less clear in cold nuclear matter experiments. 
The E325-KEK collaboration [11] has reported a measurable shift in the masses of 
vector mesons compatible with theoretical predictions based on Brown-Rho scaling and 
QCD sum rules. On the other hand, the JLab-CLAS experiment [12] has obtained results 
compatible with vanishing mass shift, as predicted by most in-medium hadronic many- 
body analysis, where broadening is the dominant effect [13, 14]. 

The c//o(600) meson is also a suitable candidate to study signals of chiral restoration 



through its in-medium behaviour, since it has the quantum numbers of the vacuum. A 
simple 0(4) model description would suggest that the dropping of (<j) ~ a/ (qq) driven 
by chiral restoration should imply a significant reduction of the o mass, at least in 
the chiral limit. Such a reduction could produce experimentally an enhancement of %% 
scattering and cross section near the point where M a — > 2m n and the phase space shrinks. 
This is the so called threshold enhancement effect [15], whose signals are observed in 
cold nuclear matter reactions %A — > tztzA' [16, 17] and yA — > %%A' [18]. 

Our approach is to consider the temperature and density modifications of the lightest 
meson resonances in the unitarized chiral lagrangian framework [19, 20, 21, 22]. We 
rely mostly on Chiral Perturbation Theory (ChPT), which provides a model-independent 
low-energy expansion for a given nn scattering partial wave with total isospin / and 
angular momentum J as t IJ (s; T) = ^(s) + t ! 4 J '(s; T) + 0(p 6 ) with s the center of mass 
energy squared and T the temperature. According to the usual ChPT scheme, ti is the 
tree level contribution from the lowest order lagrangian Jz?2, while t$ includes both 
the tree level from Jd4 and the one-loop diagrams from Jzf2- The latter contain the 
dependence with the temperature [19], as well as the imaginary part demanded by 
unitarity, which at finite temperature reads Im?4 = Orl^l 2 , where: 

C!r(j) = Ob(j)[l + 2n(v^/2)] (1) 

is the two-pion thermal phase space, with n(x) = [exp(x/T) — 1] _1 the Bose-Einstein 
distribution function and <Jq(s) = a/I — Am\js. The ChPT series cannot reproduce a 
resonant behaviour, which can be recovered by using a unitarization method, which 
amounts to find a unitary amplitude t satisfying exactly Imt = <Jr\t\ 2 and matching the 
ChPT series when expanded at low energies or temperatures. We follow the Inverse 
Amplitude Method (IAM) which reads: 

jam his) 2 

t 2 (s)-t 4 (s;T)+A(s;T) 

where A (s; T) is a known function constructed from ti and U ensuring that the zeros of 
the unitarized amplitude coincide with those of the perturbative one and which can be 
alternatively obtained using dispersion relations [23]. It is irrelevant for poles far from 
the real axis but should be included in medium analysis where, as commented above, 
chiral symmetry may induce a vanishing imaginary part and in fact the A function in 
that case prevents the appearance of spurious poles [21]. 

Once we have a unitarized nit amplitude, the o and p are generated as poles in 
the second Riemann sheet dynamically, i.e., without introducing them explicitly in the 
lagrangian nor assuming anything about their nature. The results with the finite-r IAM 
are showed in Fig. 1 . The so-called low-energy constants of ^4 are fixed to yield the mass 
and width of the p at their physical values at T = 0. The main conclusions of our analysis 
are: i) The main effect for the p is thermal broadening, in agreement with dilepton data, 
the mass being reduced only slightly for temperatures below chiral restoration, ii) such 
broadening does not only come from phase space increasing Oj / (Jo but also from the 
increasing of the effective pKK vertex, extracted from the residue of the amplitude at 
the pole position, iii) in the o channel, mass reduction dominates, presumably driven by 
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FIGURE 1. Results for the p (upper panels) and a (lower panels) thermal mass and width extracted 
from the IAM poles at finite temperature. 



chiral restoration, and at some point overcomes the thermal width phase space increasing 
so that the imaginary part (width) reaches a maximum value and then starts decreasing, 
iv) Near chiral restoration, the width of the a remains sizable (its spectral function 
is not peaked around the mass as in a Breit-Wigner resonance) so that no threshold 
enhancement is expected in a Heavy Ion environment, v) The mass dropping pattern in 
both the o and p cases does not scale with the condensate as expected from Brown-Rho 
scaling or simple 0(A) model arguments for the o, not even near the chiral limit. 

The next question is how to incorporate the effect of the nuclear density p in our 
approach. At T = 0, a simple procedure [21] which encodes most of the relevant 
corrections is to scale the pion decay constant f n according to [24] : 

fAp) (*>(P) _.( x _ ^ ( 1-0.35^ (3) 



fi(0) (qq)(0) V '"IfiWj \ ' Po 

where p is the nuclear density, o^n ~ 45 MeV is the pion-nucleon sigma term and po — 
0.17 fm -3 is the normal or saturation nuclear matter density. In this way we are actually 
examining the scaling properties of the resonances when we scale /|(p), or equivalently 
the quark condensate. Actually, in this simple approach there is no broadening source, 
so that only the masses change. The results of the mass scaling in our IAM approach 
are showed in Fig. 2. In the 7 = 7 = channel, the scaling corresponds to the a of the 
0(4) model. When the pole approaches the real axis, threshold enhancement is clearly 
observed and when it reaches it a pole-doubling occurs: one of the poles remains near 
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FIGURE 2. Scaling of the mass in the IAM, compared with that of f K {p). In the / = J = channel, the 
masses displayed correspond to the lowest masses of the second-sheet poles. 



threshold and eventually jumps into the first Riemann sheet becoming a nn bound state, 
while the other one (plotted in Fig. 2) tends to degenerate with the pion as its chiral 
partner. Although this phenomenon takes place at rather high p ~ 2po for this approach 
to be fully trusted, it has been obtained also in [25] and reveals the molecular-like 
character of the fo state [21]. The chiral degeneracy between the dynamically generated 
a and the pion for f n — >■ had been also noticed in [26]. In the p channel, the f K scaling 
induces a linear M p scaling in agreement with the experimental value of [11] and the 
theoretical predictions of [6] and [7], which however must be taken with a pinch of salt 
because all the broadening effects coming from many-body interactions are neglected. 

The picture that emerges from our analysis is then that when broadening effects are 
ignored, chiral restoration in the a channel takes place more or less along the expected 
0(4) pattern. However, the 1/0 study shows that thermal broadening distorts the 
picture of a qq-\ike narrow state and, despite experiencing a large mass reduction, the 
state is still wide near the transition. This is confirmed by the detailed many-body 
calculation performed in [21] for the o channel, including both nuclear density and 
temperature effects. The unitarization scheme used is the Bethe-Salpeter one, more 
suitable for nuclear many-body analysis, and takes into account the pion self-energy 
in the nuclear medium plus all the other relevant diagrams compatible with chiral 
symmetry at the same order in density. The results in Figure 3 for the imaginary part 
of the amplitude show some distinctive features of this full approach: the abrupt change 
with respect to the vacuum or thermal case, originated by physical in-medium excitations 
such as particle-hole or A-hole, which produce also strength below threshold, the sizable 
threshold enhancement, even though the pole is still rather away from the axis like in 
the thermal case, and finally the amplification of the threshold strength by the combined 
effect of temperature and density, which in fact accelerates the migration of the pole 
towards threshold by sudden mass decrease. 

PIONS OUT OF CHEMICAL EQUILIBRIUM 

In the hadron gas formed in a Relativistic Heavy Ion Collision, it is phenomenologically 
well justified to consider a phase where thermal equilibrium still prevails but inelastic 




FIGURE 3. Imaginary part of the kk amplitude in the I = J = channel at finite temperature and 
nuclear density in the many-body BS approach. 



collisions involving particle number change are negligible, so that chemical equilibrium 
is lost. For instance, according to the estimates in [27], at T = 150 MeV the relaxation 
time of elastic nil collision is X e i ~ 2fm/c, whereas that of the process %% <-> %%%% is 
Tin ~ 200 fm/c, the hadronic phase lifetime being about 10 fm/c. In other words, there is 
a temperature window between the chemical freeze-out temperature T cnem ~ 180 MeV, 
below which approximately only elastic collisions remain, and the thermal freeze-out 
one Tfher ~ 100-120 MeV of hadron decoupling, where the chemical potential associated 
to pion number conservation \i n is not zero [28, 29, 30]. Actually, in this phase \x % is a 
function of T so that — > T c hem) = 0. On the other hand, different phenomenological 
fits of particle yields and ratios at SPS and RHIC energies predict \ijiijtha) ~ 70-100 
MeV [29, 30]. 

In a recent work [31], we have considered a ChPT description of this phase in the pion 
gas, within a Quantum Field Theory treatment of chemical nonequilibrium for boson 
fields. A crucial point is that the total pion number is only approximately conserved, 
unlike for instance the electric charge or the third isospin component. This makes the 
usual path-integral description in terms of field states inconvenient for this case. We have 
used instead an holomorphic representation, which allows to develop the formalism in 
terms of creation and annihilation operators, in terms of which one can easily express 
the free particle number. Another technical complication is that there are no KMS-like 
periodicity conditions, which prevents the imaginary-time thermal formalism to be used. 
As in other nonequilibrium formulations, the appropriate choice is a real-time contour 
in complex time. This implies the doubling of the pion degrees of freedom, so that the 
propagators have a two-dimensional matrix structure. One of the two types of fields 
is unphysical and only appears in internal lines. After a detailed analysis, we obtain the 
generating functional, which can be read from the usual one at \l n = with the following 
replacements for the free propagators and for the free partition function: 
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with V the system volume and E 2 = \p\ 2 + With the above ingredients, we can 
calculate the different thermodynamical variables of the pion gas in this regime, by 
evaluating the corresponding closed diagrams, whose topology is the same as the \l n = 
case [32]. We have extended our analysis up to £?(T 8 ). At that order, particle-changing 
processes can be seen to become already important, but the perturbative scheme remains 
consistent since an additional physical condition has to be imposed in order to regulate 
the plasma dynamical evolution from thermal to chemical equilibrium, namely the 
/x^(r) function. A simple way to model this is to demand that the ratio of entropy 
density to pion particle number density s/n remains constant [28]. This isentropic 
condition is nothing but the combination of entropy conservation without dissipation 
and approximate particle number conservation in this phase. The crucial point is that 
this ratio is a decreasing function of both T and \i n , so that as the systems cools down, 
the chemical potential grows in order to keep it constant. The precise value of s/n can 
be fixed either by the chemical freeze-out temperature at which \l%(J c hem) = or by the 
value \x n (J t h er ) at thermal freeze-out. We show in Figure 4 our results for different chiral 
orders, comparing also with the virial expansion approach followed in [33]. We see that 
one of the consequences of including the pion interactions is to lower the value of T c i iem 
by about 25 MeV with respect to the ideal gas. 

Our diagrammatic approach allows also to calculate the pion self-energy corrections 
in T and The detailed analysis and results are given in [31]. To leading order in the 
pion density (dilute gas regime) the self-energy obeys a Luscher-like relation: 
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where tI k (s) is the forward kk scattering amplitude. The 0(p 4 ) and unitarized (IAM) 
amplitudes produce a decreasing pion mass with both T and \x n . Actually, at a given tem- 
perature, the mass can reach the chemical potential, as showed in Figure 4. This implies 
the interesting possibility that pion Bose-Einstein condensation could be reached dy- 
namically, i.e., driven by interactions. The corresponding /i BE (T) where this condition 
is met is above but not far from the isentropic curves for realistic chemical freeze-out 
conditions. On the other hand, from the pion width T p , we can estimate the elastic mean 
collision time x e \ and the thermal freeze-out temperature where T e /(7^ er ) ~ 10 fm/c. 
Taking into account the \l n dependence of T e / through the isentropic yields a 

reduction KT ther ~ -20 MeV with respect to neglecting chemical nonequilibrium. 



Transport coefficients measure the linear response of the system to a deviation from 
equilibrium induced by an external source. Their knowledge is essential to describe 
many phenomenological aspects regarding dissipation effects in Heavy Ion Collisions. 
Recently, we have developed an extensive programme for the calculation of transport 
coefficients within diagrammatic ChPT, with direct application to the mesonic phase of 
the plasma expansion [34, 35, 36]. The diagrammatic treatment is technically nontrivial 
because transport coefficients involve the zero external frequency and momentum limit 
of retarded correlators. This implies the appearance of the so called pinching poles, 
proportional to 1 /T with T the collisional width of the particles in the medium. This 
generates nonperturbative contributions, since T is proportional to the collision ampli- 
tude (see e.g. eq.(5)) and therefore perturbatively small. Thus, the usual chiral power 
counting of ChPT has to modified to account properly for these contributions [34] . 

The viscosity coefficients are particularly interesting for phenomenology. The shear 
viscosity to entropy density ratio rj/s can be determined through elliptic flow analysis, 
the present RHIC data pointing to a small value r\/s < 0.5 [37], an almost perfect 
fluid. Our ChPT result is greatly influenced by a correct unitarized description of the 
scattering amplitudes entering the pion width. In fact, without unitarization, at <f?(p 4 ), 
rj/s is a monotonically decreasing function of T, violating the so called KSS bound 
rj/s > \/Ak [38], as showed in Figure 5. The unitarized result satisfies that bound, gives 
phenomenologically compatible values at the relevant temperatures and is compatible 
with the existence of a minimum for this ratio near the transition [39]. The case of 
the bulk viscosity £ is also important because it is directly sensitive to variations of 
conformal symmetry. Thus, it should have a similar behaviour [40] to the trace anomaly: 



where Ty is the energy-momentum tensor and P the thermodynamic pressure. The trace 
anomaly has a maximum at the QCD transition, originated mostly from the anoma- 
lous glue contribution to the beta function and observed in the lattice [2]. Therefore, 
significantly large values for £ could be observed near the critical point. Our results 
in Figure 5 show a clear correlation between these quantities, both developing a two- 
peak structure. The first maximum comes from the explicit conformal breaking due to 
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the quark mass and thus it disappears in the chiral limit. Actually, in that temperature 
range 8(0)t ~ 2m q d{qq)j [35]. The second maximum is the critical one and is al- 
most unchanged in the chiral limit, not being related to the quark condensate or to the 
quark mass, but to deconfinement effects encoded in the gluon condensate anomalous 
conformal breaking. We see again that incorporating properly unitarization is crucial to 
describe correctly the high temperature region. It is also essential to consider the speed 
of sound squared cl to <^(r 8 ) in order to capture its minimum at the transition aris- 
ing from anomalous conformal breaking. Our analysis does not rely on assumptions on 
the behaviour of the 60 spectral function, as it is the case of the original proposal in 
[40], which has been recently corrected [41]. It is therefore unclear from general con- 
siderations how the correlation between the trace anomaly and the bulk viscosity should 
actually hold. In fact, it does not hold in simple models mimicking QCD [42]. We be- 
lieve then that our pion gas description can be useful to understand this and other issues 
regarding the behaviour of transport coefficients below the critical temperature. 
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